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Abstract 

The electronic structure and magnetic properties of CrSb2 have been investigated by ab-initio 
calculations with an emphasis on the role of the magnetic structure for the ground state. The 
influence of correlation effects has been investigated by performing fixed spin moment (FSM) 
calculations showing their important role for the electronic and magnetic properties. The details 
of the electronic structure of CrSb2 are analyzed by a comparison with those of FeSb2 . The results 
obtained contribute in particular to the understanding of the temperature dependence of transport 
and magnetic behavior observed experimentally. 

PACS numbers: Valid PACS appear here 
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I. INTRODUCTION 



According to the experimental data available in the literature, the CrSb phase diagram 
shows two stable phases, having NiAs (CrSb) and marcasite (CrSb2) structures lM3|, re- 
spectively. In the case of CrSb2 with marcasite structure (space crystal group Pnnm) the Cr 
atoms are octahedrally coordinated by six nearest neighbor Sb atoms. As discussed in the 
literature, e.g. 4J, CrSb 2 belongs together with FeSb 2 to the so-called 'Jahn- Teller' or 'regu- 
lar' class of marcasites (characterized by a ratio of the lattice parameters c/a ~ 0.53 — 0.57, 
c/b « 0.48) in contrast to 'anomalous' marcasites (characterized by c/a ~ 0.73 — 0.75, 
c/b ~ 0.62) [4]. The 'regular' marcasites have an electronic structure leading to a Jahn- Teller 
instability with respect to octahedral symmetry, i.e. to a compression of the octahedral sur- 
rounding. This leads to a C^h symmetry at the transition metal (TM) sites that gives rise 
to a corresponding splitting of the d-states in the crystal field {4, 5|. According to the ele- 
ment specific occupation of the splitted d-shell one can distinguish between the properties 
of low-spin 3d 4 configuration occurring for FeSb2 and high-spin 3d 2 configuration of CrSb2- 
Although the CrSb2 and FeSb2 compounds differ concerning these features, they have very 
similar crystal structures and therefore, it seems worth to consider them in parallel. In 
particular, is also motivated by the fact that only little information on CrSb2 is available 
while the properties of FeSb2 are investigated experimentally and theoretically in an ex- 
tensive way [q-lllj. From the literature it is known that FeSb2 and CrSb2 compounds are 
narrow-gap semiconductors at ambient pressure [l 3]]. While FeSb2 is nonmagnetic (NM) 
at low temperature, CrSb 2 has an antiferromagnetic (AFM) order with a magnetic moment 
of the Cr atom of about 1.94/ig [12I. \l3\. A temperature increase results in a transition from 
the small gap semiconductor to a metallic state with strong spin moment fluctuations {(J. 
This behavior is attributed to strong electron-electron correlation effects [(3, 0-1(3, [ijj] that 
gives rise to the properties observed experimentally: the temperature dependence of mag- 



netic susceptibility 



the temperature dependent behavior of the electrical resistivity p(T) 
15l . Il6| and colossal thermopower S(T) at T = 10 K [15|, ll6(. This behavior was found to 
be similar to that observed in FeSi with a Kondo insulator (KI) model as a possible scenario 
to describe the observed properties which have been investigated theoretically via LDA+U 
electronic structure calculations Jl7j . The corresponding investigations on FeSb2 have been 
performed by Lukoyanov et al. [9j giving arguments for the KI description of the spin state 
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in FeSb2- However, the origin of the electronic properties of FeSb2 as well as of FeSi, is still 
under discussion. In particular the recent theoretical investigations [h], [ijj gave arguments 
against the KI scenario of transition in these systems and discuss the effect of strong local 
dynamical correlations. 

As was mentioned above, CrSb2 has an AFM order below T 273 K with a local magnetic 
moment of Cr atom of about 1.94//B [12I Q]- The heat capacity observed in experiment has 



a lambda-shaped maximum at T/v = 274.08 K [18(, which corresponds to a transition to the 
AFM state upon temperature decrease. The magnetic susceptibility x(T) has an anomaly 



at T = 55 K which may be attributed 



19] to a crossover between the different electronic 



states in the AFM semiconductor, but cannot be related to any magnetic transformation. 



No anomaly of the magnetic susceptibility has been observed near T/v = 273 K [19] . On the 
other hand, the electrical resistivity measurements show an anomaly at the Neel temperature 
273 K, and at the same time has a plateau in the temperature regime 50 — 80 K. This finding 
is supposed to have a common origin with the susceptibility anomaly and is also attributed 
to a transition between the different semiconductive electronic states with the energy gap 
evaluated as A = 0.07 eV. The behavior sketched for the resistivity has been observed by 



various experimental groups 
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19l-|21|. In addition, Hu et al. 13J have demonstrated a 



rather pronounced anisotropy in the transport properties of FeSb2 which is not observed in 



CrSb 2 . The measurements of the transverse magnetoresistance (TM) 



p(h)-p(Q) 
ftp) 



in a magnetic 



field up to H = 90 kOe at different temperatures exhibits an increase of this value with H 



and then linearly as H increases further 



19]. As has been pointed out [191 ]. a similar behavior 



of the TM has been observed also for the FeSi semiconductor. The temperature dependent 
behavior of exhibit a maximum at around ~ 20K. 



Li et al. 



21] have shown that the thermal conductivity of CrSb2 has a maximum at 



T w 50 K, that they attributed to the effect of phonon-phonon Umklapp scattering. The 
thermopower S measured as a function of temperature is negative in the whole temperature 
range and has a large peak at ~ 60 K which correlates with the position of the plateau of 
the temperature dependent resistivity p(T). For comparison, the thermopower \S\ in FeSb2 



m 



160 



has a giant maximum at T « 10f^ and its thermal conductivity at « 15 K. 

A high-pressure experiment on CrSb2 showed that for a pressure above 5.5 GPa a phase 
transition to the CuAl 2 structure occurs accompanied by ferromagnetic (FM) order with a 
Curie temperature T « 160 K and a saturation magnetic moment of about 1.2/xb per Cr 
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atom [22j. Experiments on FeSb2 show no phase transitions up to 7 GPa [6]. This is in line 
with theoretical investigations that predict that for a transition from the marcasite to the 
CuAl 2 phase the pressure of about 41 GPa is required [23 ]. 



II. COMPUTATIONAL DETAILS 



The calculation are done with the ELK code based on the FP-LAPW method [24J. To 
take into account correlation effects for Fe and Cr the calculations were done within the GGA 
as well as the GGA+U method. We used the exchange correlation functional by Perdew, 
Burke, Ernzerhof (PBE) (2^]. For the GGA+U calculations the Coulomb (U) and exchange 
(J) parameters were chosen as follows: U = 2.7 eV and J = 0.3 eV for Cr and U = 2.6 eV 

n 

and J = 0.88 eV for Fe [9J. As double counting correction the around mean field expression 
(AMF)|26-28] was used. The Cr magnetic moments presented below are calculated within 



the MT sphere with radius Rmt = 2.32 a.u. The experimental structure parameters used 
for CrSb2 with marcasite structure are: a = 6.1481 Lb/a = 1.1404, c/a = 0.5428 [29(, and 
for FeSb 2 : a = 5.8328 A, b/a = 1.1208, c/a = 0.5482 [29j. 



III. RESULTS 



Although there are several experimental results in the literature on CrSb2, no detailed 
theoretical investigation have been done so far. Despite clear differences with FeSb2, the 
two systems have also some common features. Therefore, it is interesting to analyze the 
properties of CrSb 2 in comparison with those of FeSb 2 . 

As a first step, electronic structure calculations for FeSb2 and CrSb2 have been performed 
based on the GGA using the experimental structure parameters. In the case of CrSb2 with 
a small unit cell these calculations lead to a FM metallic ground state of the system (see the 
corresponding density of states (DOS) plot on the Cr atoms in Fig. [lb), in contradiction 
to the experimental results that show a small-gap semiconducting behavior for the system 
with AFM order. To describe the AFM order observed experimentally, supercell calculations 
are required, which will be discussed below. As a first step, we focus on the FM state of 
the system. In comparison, GGA calculations for FeSb2 lead to a semiconducting NM 
state with the Fe DOS shown in Fig. [T^,. GGA+U calculations performed to explain the 
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FIG. 1: Density of states on the transition- metal atoms in FeSb2 (NM state) (a), and in CrSb2 
in NM (GGA results) and weakly magnetic with constrained magnetic moment mc r = 0.05^# 
(GGA+U results) (b), FM (c) and AFM (d) states. Solid and thin dashed lines represent the 
results of GGA+U (U = 2.7 eV for Cr and 2.6 eV for Fe) and GGA calculations, respectively. 
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FIG. 2: Electronic band structure for NM FeSb 2 via GGA (a) and for FM CrSb 2 via GGA (b) and 
GGA+U, U = 2.7 eV (c) 



finite temperature properties of FeSb2 as reported in the literature demonstrate the crucial 
importance of a proper description of local correlations this compound. Motivated by this 
experience, GGA+U calculations have been performed also for CrSb 2 . 

To investigate the dependence of the results on the U parameter the calculations have 
been performed for different U value varying from eV to 7 eV (see, e.g., the discussions 



m 



171). In the following, if not otherwise noted, the results are presented for U = 2.7 eV as 



this value gives a well defined semiconducting state for CrSb 2 with FM order. However, it 
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should be noted, that in the case of AFM order, that is found experimentally jl^, Q for to 
the ground state, the U value does not play a crucial role for the semiconducting properties. 

Fig. [Tb displays the DOS on the Cr sites for NM CrSb2 exhibiting metallic behaviour for 
GGA calculations (dashed line) with an energy gap above the Fermi level. Local correlations 
accounted for via GGA+U approach (with U = 2.7 eV) lead to a pronounced modification 
of the electronic structure even for a slightly spin-polarized system within LSM calculations 
with Cr magnetic moment of 0.05/ib (Fig. [lb, solid line). This results in a shift of the Fermi 
energy into the energy gap, i.e. the system becomes semiconducting. 

In the case of CrSb 2 with FM order, selfconsistent GGA calculations (see the DOS in 
Fig. [Tb) lead to a Cr magnetic moment of about 2.63/i# and a metallic state of the system. 
GGA+U results on the other hand exhibit the energy gap at the Fermi energy, different 
for the two different spin channels, and a magnetic moment of Cr atom of about 1.64//g. 
Corresponding modifications of the electronic structure can be seen in the energy band 
dispersion shown in Fig. EJ b (GGA) and c (GGA+U). 

Calculations for the AFM state of CrSb2 have been performed for the magnetic structure 
obtained experimentally [l^j]. Interestingly, in this case a small energy gap at the Fermi 
level occurs even within the GGA calculations. Local correlations treated via GGA+U only 
enhance the energy gap (see Fig. [TH) but do not result in further noteworthy modifications 
of the electronic structure. 

Total energy calculations as a function of the volume have been performed for the NM, 
FM and AFM states of CrSb2 with marcasite structure as well as for CrSb2 in the FM state 
with CuAl2 structure. The results are presented in Fig. [U The volume in these calculations 
have been varied by a variation of the lattice parameter a, keeping the ratios b/a and c/a 
fixed as taken from experiment: b/a = 1.1404 and c/a = 0.5428 for marcasite structure 29] . 
and c/a = 0.885 for CuAl 2 structure 22J. All calculations have been performed within the 
GGA+U (U = 2.7 eV) approach. Note that in the case of the CuAl2 structure the difference 
in DOS results obtained via GGA and GGA+U is rather small and in both cases a metallic 
state of the system is observed (see Fig. [3]). 

From these calculations one can see that the AFM state for CrSb 2 with marcasite struc- 
ture is indeed the ground state of the system as found in experiment, while the FM and 
NM states have slightly higher total energies. All total energy curves have a minimum at 
the volume which is nearly the same for different types of magnetic order and fits well to 
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FIG. 3: Density of states on the Cr atoms for CrSb2 with CuAi2 structure in FM state. Solid and 
dashed lines represent the results of GGA+U (U = 2.7 eV) and GGA calculations, respectively. 
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FIG. 4: Total energy of CrSb2 calculated using the GGA+U for various magnetic phases with mar- 
casite structure and for FM phase with CuAi2 structure. The c/a and c/b ratios for all calculations 



have been kept constant at the experimental value 



22 



the experimental value. The curves for the marcasite structure have a crossing point with 
the curve for the CuAl 2 structure that has a minimum at a smaller volume than in the case 
of marcasite structure. These results are fully in line with the experimental data showing a 
transition from the marcasite to the FM CuAl2 phase at the pressure of about 5.5 GPa 22 ]. 

The temperature dependence of transport and magnetic properties of CrSb2 and the 
influence of correlation effects can be analyzed by performing fixed spin moment (FSM) 
calculations and comparing the results with the results for FeSb2- As was reported, e.g., by 
Lukoyanov et al. [9|, the ground state of FeSb 2 is NM (see Fig. [5^). GGA+U calculations, on 
the other hand, result in an additional minimum at (Fig. 0d) that is used to explain the 
strong spin fluctuations at high temperatures. According to the discussions by Lukoyanov 
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FIG. 5: The results of FSM calculations for FeSt)2: Total energy per unit cell vs fixed Fe magnetic 
moment via GGA (a) and GGA+U (U = 2.6 eV) (b). 
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FIG. 6: Fe DOS in FeSb 2 obtained via GGA+U: for the NM (a, U = eV) and FM (b, U = 2.6 
eV) states. Thick solid line corresponds to 3d-states with d z 2 character, the dashed line shows the 
DOS of states with d x 2_ y 2 character. 

et al. [9] on FeSb2 as well as Anisimov et al. [17] on FeSi, this minimum is associated with 
the occupation of narrow Fe 3d-bands with d z 2 character at the bottom of the conduction 
band, and states of d x 2_ y 2 character near the top of the valence band (see Fig. E]). The 
results of corresponding FSM calculations for CrSb 2 are shown in Fig. [7J They exhibit 
rather different features when compared to FeSb2. The GGA calculations show an absolute 
minimum for the total energy around a Cr magnetic moment of 2.0/iB (Fig- EK), in contrast 
to FeSb2 exhibiting a NM ground state. This means, the NM state of CrSb2 (with the Cr 
DOS shown in Fig. [Ht) is unstable with respect to the creation of a local magnetic moment 
on the Cr atoms. In general, this can lead either to FM or AFM order in the system with the 
corresponding Cr DOS shown in Fig. [HI b and d, respectively. This minimum is split into a 



deeper at 1.8/xb and a more shallow one at ~ 2.6//£ which depends essentially on the U 
value. The Cr DOS curves corresponding to these two minima are plotted in Fig. [9] a and 
b. One can see that the state with 1.8/i# is associated with the unoccupied majority-spin 
Cr 3<i 2 2-states, while magnetic moment « 2.6/xb is caused by additional occupation of Cr 
states of d z 2 character. 

In the case of a higher magnetic moment the exchange splitting of the electronic states 
leads to an occupation of antibonding states in the spin-up channel and depleting of the 
bonding states in the spin-down channel. Such an occupation leads to a partial compensation 
of the energy gain due to the exchange splitting. Therefore, the FM state prefers to have an 
unoccupied 3d z 2 spin-up energy band, that is stabilized by treating local correlation effects 
via the GGA+U approach (see the DOS for U = 2 eV in Fig. c and d). Thus, in the case 
of U = 2eV the total energy within the FSM calculations has one minimum at m = O/ig 
and a second broad minimum in the regime 1.7 to 2.0/ie- Fig. [9] c and d shows the DOS 
corresponding to the values of Cr magnetic moment m = 1.7 {is and m = 2.0fis- These 
results (for U = 2 eV, 2.7 eV) exhibit a similar behavior as in the case of FeSb2 with the 
difference that the minimum at the finite magnetic moment is deeper. A further increase of 
the U value above 5. eV results in an instability of the FM state. 

Another scenario is possible with a shift downwards in energy of the majority-spin 3d z 2- 
states of Fe, changing their relative position with respect to the the p states of Sb. This also 
leads to the arrangement of the Fermi level within the energy gap. This scenario is associated 
with an AFM state with the corresponding DOS shown in Fig. [HE (for U = 0) and Fig. 
Efe (for U = 2.7 eV). In this case GGA calculations result in a small-gap semiconducting 
state with a local Cr magnetic moment close to that obtained for the FM state (see Table 
1), since Cr spin- up 3<i z 2-states in this case are also occupied. In contrast to FM case, the 
energy gap at Ep in this case occurs due to different spin-dependent hybridisation of the Cr 
d states with p states of Sb atoms. Also, in contrast to the FM state treated via GGA, AFM 
order in CrSb2 is expected to be more stable since all bonding electronic states are occupied 
and antibonding states are not occupied. Altogether the GGA+U calculations result only 
in a slight variation of the electronic structure leading in particular to a decrease of the Cr 
magnetic moments. 

From the results above one can draw the following conclusions. Stabilization of the AFM 
order as a ground state of CrSb2 is related to a shift below the Fermi energy of Fe spin-up 
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FIG. 7: The results of FSM calculations for CrSb2 with FM order. Total energy per unit cell vs 
fixed Cr magnetic moment via GGA (a) and GGA+U: U = 2 eV (b), U = 2.7 eV (c),U = 5.0 eV 
(d). 
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FIG. 8: Cr DOS for CrSb 2 in NM state (a, U = 0.0 eV), FM state (b, U = 0.0 eV) and (c, U = 2.7 
eV), AFM state (d, U = 0.0 eV) and (e, U = 2.7 eV). Thick solid line corresponds to the 3d-states 
with d r 2 character. 



11 



rP-l 



■3 0.5 



a -0.5 
fa 




-2 
Energy (eV) 



-2 
Energy (eV) 



U 0.5 

> 

"3 
55 

X -0.5 
w 

d -1 



(a) 



1 1 1 1 
" m=1.7Li (U=2.0eV) . 


1 A 1 1 


i,i, 




1 , 1 



£ 0.5 

> 

a 

& -0.5 



(c) 



1 1 1 1 
- m=2.6|i B (U=0 eV) 










: 




i 



-2 
Energy (eV) 



-2 
Energy (eV) 



(b) 



11 

m=2.0ji B (U=2.0 eV) 

— 


A ■ 




i 



(d) 



FIG. 9: Results of FSM calculations for CrSb2- Cr DOS in the FM state for the constrained Cr spin 
moment m = ( a ) an d m = 2.6[Ib (b) for U = 0.0 eV); and m = 1.7//^ (a) and m = 2.0/i^ 

(b) for U = 2.0 eV). Thick solid line corresponds to the 3d-states with d z i character. 

states of d z 2 character with a corresponding modification of the spin-dependent hybridization 
of Cr d states and p states of Sb. This moves the Fermi level into the energy gap to make 
all bonding states to be occupied and antibonding states unoccupied. Local correlations 
treated via GGA+U influence in this case mainly the width of the energy gap. In the case 
of the metastable FM state GGA+U calculations lead to more crucial modifications of the 
electronic structure resulting in a semiconducting state of the system. This occurs due to 
a push up of the majority-spin d z i states of Cr above Ep in cost of stronger hybridization 
of minority-spin states with p states of Sb keeping the total number of electrons on the Fe 
atoms unchanged. 

Although the FSM calculations have been performed only for the FM order, the results 
above allow also to expect that being in the AFM state, CrSb2 should exhibit a metastable 
state determined by a different occupation of the d z 2 Cr states. Indeed, a temperature 
induced transition to the metastable state could explain the behavior of the susceptibility 
and transport properties observed experimentally at low temperature. It is possible that 
these two AFM states with slightly different electronic structure can be characterized also 
by different Cr local magnetic moments. 
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IV. CONCLUSION 

In summary, one has to stress that the AFM order in CrSb 2 plays a crucial role for the type 
of conductivity in the system exhibiting the small-gap semiconducting properties. Although 
the local correlations treated via the GGA+U approach influence mainly the width of the 
energy gap of CrSb2, their role can be important to describe the experimental results as it 
was shown for FeSb2- In general, the present theoretical results are in a good agreement 
with the experiment, showing the AFM ground state of CrSb 2 with marcasite structure, 
with a Cr local magnetic moment of about 2.03/zb and with a small energy gap at the Fermi 
level. Concluding from the total energy calculations a phase transition under pressure to 
the CuAl 2 phase exhibiting the properties of FM metal should be observed. 
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